Vortex Interactions and Thermally Induced Crossover from Type-I to Type-II 

Superconductivity 

J. Hove,! S. Moj] and A. Sudb0| 

Department of Physics 
Norwegian University of Science and Technology, 
N-7491 Trondhcim, Norway 
(Dated: February 1, 2008) 

We have computed the effective interaction between vortices in the Ginzburg-Landau model from 
large-scale Monte-Carlo simulations, taking thermal fluctuations of matter fields and gauge fields 
fully into account close to the critical temperature. We find a change, in the form of what appears 
to be a crossover, from an attractive to a repulsive effective vortex interactions in an intermediate 
range of Ginzburg-Landau parameters k G [0.76, 1]/\/2 upon increasing the temperature in the 
superconducting state. This corresponds to a thermally induced crossover from type-I to type-II 
superconductivity around a temperature Tc-r{i^), which we map out in the vicinity of the metal- 
to-superconductor transition. In order to see this crossover, it is essential to include amplitude 
fiuctuations of the matter field, in addition to phase-fluctuations and gauge-fleld fluctuations. We 
present a simple physical picture of the crossover, and relate it to observations in Ta and Nb 
elemental superconductors which have low-temperature values of n in the relevant range. 
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I. INTRODUCTION 

The nature of the phase-transition in systems of a 
scalar matter field coupled to a massless gauge-field has 
a long history in condensed matter physics, dating at 
least back to the introduction of Xhe Ginzburg-Landau 
(GL) theory of superconductivityfiJ. At the mean-field 
level, ignoring spatial variations in gauge fields as well 
as matter fields leads to the prediction of a second order 
phase transition in the model, with classical mean-field 
exponents for all values of the GL parameter k. The 
first attempt to seriously consider the role of fluctuations 
on the order of the metal-to-superconductor transition 
was made by Halperin, Lubensky, and MacI, who found 
that ignoring matter-field fluctuations entirely, and treat- 
ing gauge-field fluctuations exactly, resulted in a perma- 
nent first order transition for all values of k, since the 
gauge-field-fluctuations produced an extra term ~ — 10|^ 
in the matter field sector of the theory in three spatial 
dimensions, where the complex matter field is denoted 
by (j) and represents the condensate order parameter. (In. 
the context of particle physics, Coleman and Weinbergcl 
studied the equivalent problem of spontaneous symme- 
try breaking due to radiative corrections in the Abelian 
Higgs model in four space-time dimensions, finding the 
additional term 0"* ln(0^/(/)Q), where the real rriatter field 
is denoted 4> and represents a scalar meson. )□. Subse- 
quently, Dasgupta and HalperinEl found, using duality 
arguments in conjunction with Monte-Carlo simulations, 
that when gauge- field fluctuations and phase fluctuations 



of the scalar matter field are taken into account, but am- 
plitude fluctuations are ignored, the phasje transition is 
permanently second ordertl. BartholomewQ then reported 
results from Monte-Carlo simulations for the case when 
also amplitude fluctuations are taken into account, con- 
cluding that the phase transition changes from first to 
second order at a particular value of the GL parame- 
ter K « As far as this numerical value is con- 
cerned, note that the problem of finding a tricritical 
value Ktri separating first and second order transitions 
is extremely derrLanding even by present day supercom- 
puting standardsQ (jsee below). Using ingenious duality 
arguments, Kleinerta obtained that the change from first 
to second order transition should occur at k w 0.8/\/2. 
The value of k that separates a first order (discontinuous) 
transition from a|Second order (continuous) one, defines 
a tricritical pointlj, and will hereafter be denoted Ktri- 
Note that to obtain the above result, it is necessary to 
allow for amplitude fluctuations in the superconducting 
order parameter, which become important for small to 
intermediate values of k, but are totally negligible in the 
extreme type-II regime k >> 1. 

The critical properties of a superconductor may be 
investigated at the phenomenological level by the GL 
model of a complex scalar matter field 4> coupled to a 
fiuctuating massless gauge field A. It is this feature of 
the gauge field that makes the GL model so difficult to 
access by the standard techniques employing the renor- 
malisation groupaO. The GL model in d spatial dimen- 
sions is defined by the functional integral 
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where F^^ = S^Aj, — d^A^, q is the charge couphng the 
condensate matter field cj) to the fluctuating gauge field 
Afj_, A is a self coupling, and m'^ is a mass parameter 
which changes sign at the mean field critical temperature. 
When all dimensionful quantities are expressed in powers 
of the scale represented by q^, the GL model may be 
formulated in terms of the two diniensionless parameters 
y = /q^ and x = \/q^. In this case, y is temperature 
like and drives the system through a phase transition, and 
a; = is the square of the Ginzburg-Landau parameter. 
Depending on the value of x, the transition ispcither first 
order for x < Xtri, Qi continuous for x > XtriElCHS. 

In a recent paperQ, we have determined Xtri — 0.295 ± 
0.025. This corresponds to a tricritical value of the GL 
parameter Ktd = (0.76 ± 0.04)/\/2, in rather remark- 
able agreement with the results of Ref. |[ Moreover in 
Ref. it was also argued that this value of x or k is the 
demarkation value which separates type-I and type-II su- 
perconductivity, rather than the classical mean-field value 
K — I/V2. The connection can be made when one real- 
izes that criticality at the metal-to-superconductor tran- 
sition requires that topological defects of the matter field 
in the form of vortex loops are stable. On the other hand, 
there is a connection between critical exponents a»d geo- 
metric properties of a tangle of such vortex loopstil. The 
fractal dimension Dh of the vortex- loop tangle-, is con- 
nected to the anomalous scaling dimension r^^llS of the 
matter field in a field theory of ther-wpttex-loop gas, a 
theory dual to the original GL theoryt3E3 by the relation 
Dh + Vip = 2. Since the anomalous scaling dimension is 
connected to the order parameter exponent /3 of t>ip dual 
matter field by the relation 2/3 = - 2 r/^)^^, it 
follows that a collapse of the vortex-loop tangle implies 
Dh = d and hence /3 = indicative of a first order tran- 
sition. Here, d is the spatial dimension of the system. 
Now, a collapse of the tangle in turn implies an effec- 
tive attraction between vortices, or type-I behavior. On 
the other hand, a stable vortex-loop tangle at the critical 
point, with fractal dimension Dh < d, implies first of all 
type-II behavior, but also rj^ > 2 — d and (3 > 0, and 
hence a second order transition. 

The above assertion, that the tricritical value of 
K separates first order from second order metal-to- 
superconductor transition, and moreover also separates 
type-I from type-II behavior when the system is on the 
phase-transition line yc{x), is in contrast to the conven- 
tional wisdom that type-I and type-II superconductiv- 
ity is separated hy x = 0.5. Based on the above argu- 
ments, we have proposed the phase diagram shown in 
Fig. 1 of Ref. which contains a new line separating 
type-I and type-II superconductivity. The shape of this 
line was inferred from the observation that far from the 



phase transition, mean-field estimates of the boundary 
between type-I and type-II should be precise, and hence 
this boundary should asymptotically approach a; = 0.5 
from below as the temperature is reduced. 

It is the purpose of this paper to show directly, by com- 
puting the effective thermally renormalized interaction 
between vortices via large-scale Monte-Carlo simulations, 
that this quantity changes from being repulsive to attrac- 
tive in the intermediate regime k € [0.76, l]/-\/2. Since 
the sign of the vortex-interaction is the microscopic diag- 
nostics, in terms of vortex degrees of freedom, for distin- 
guishing type-I from type-II superconductivity, the large- 
scale simulations we present in this paper confirm the 
above conjectures and plausibility arguments of Ref. 0. 

In an external field the GL model has, classical so- 
lutions in terrns_of Abrikosov flux tubesEj, or Nielsen- 
Olesen vorticesEj, and the concept of type-I versus type- 
II superconductivity is based on the interaction between 
these vortices. For type-I superconductors they attract 
each other, whereas for type-II sueerconductors the in- 
teraction is repulsive. Abrikosovlla showed that at the 
mean field level type-I and type-II superconductors are 
separated at k = l/v^. We will refer to the value of n 
separating type-I from type-II behavior at kj/h, which 
we find varies with y. It is not a sharply defined quan- 
tity, since it represents a crossover line. The excep- 
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tion IS at y = 

ratp-yalculations of vortex interactions have been carried 
outlijME3ElE2l, but none of these approaches take ther- 
mal fluctuations into account. A recent overview of su- 
perconductors with K close to 1/ can be found in Ref. 
p^ , see also Ref. |2^ . 

Superconductors with k kXJ\/2 were studied exten- 
sively in the 1960s and 1970gE3, and in particular mea- 
surements on the metals Ta and Nb demonstrated that 
the noticjU-of a temperature independent value of was 
incorrectcJ. At the time, this was explained foAh a mean- 
field theory involving three GL parameter£j. Thermal 
fluctuations, not addressed at the mean-field level, offer 
an alternative and above all simpler explanation for the 
observations of crossovers from type-I to type-II behav- 
ior in one and the same compound as the temperature is 
increased. 

We have performed large scale Monte Carlo (MC) sim- 
ulations on the lattice version of Eq. |l|, with two vortices 
penetrating the sample in the z direction. By measuring 
the interaction between these two vortices we have de- 
termined the value of kj/u, in particular how this value 
is affected by thermal fluctuations close to the critical 
point. 
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II. MODEL, SIMULATIONS AND RESULTS 

To perform simulations on Eq. 0, we have defined a discrete version as foUowJl^: 



x,i<j 



In Eq. H Ui (x) — aqAi (x) and aij — ai (x) + aj (x + i) — aj (x) — aj (x + j) . Pg and [32 are related to the contimmm 
parameters x and y and the lattice constant a, 



/3g 



/32 = -^ 

PG 



y 3.1759115(1 + 2a;) (-4 + 8a; - Sx^) (In6/3G + 0.09) - 1.1 + 4.6a; 
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Note that (32 contains the effect of ultraviolet renormal- 
ization_iii.the continuum limit when the lattice constant 
a OQEEI. The model Eq. || is defined on a numerical 
grid of size x Ny x Nz, corresponding to a physical 
size of Lj. X Ly X L2, with Li ~ Nia. All our simulations 
have been on cubic systems with Pg —i^ 

To impose an external magnetic fieldEj, we modify the 



J 



action Eq. ^ by changing the field energy along one stack 
of plaquettes located at , j/g in such a way that the 
action is minimized for ai2(a;o, yo, z) = —2Trn instead of 
zero, corresponding to forcing a number of n flux quanta 
through the system. Hence, the action [a, i/"; for n 
flux-quanta forced through the system, is given by 



S [a,ip; n] ^ S [a,ip; 0] + E (27mcki2(a;o, yo, z) + 27r 
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The second term in Eq. |^ corresponds to forcing a flux 
<f>_B through the lattice in the negative z-direction 
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The crucial point is that, due to periodicity, the total 
flux through the system must be zero, i.e. 



E"i2(a;,y, z) = Vz. 



x,y 

Consequently, the n flux-quanta of the total flux 21111 /q 
must return in the +z direction. This flux returns in a 
manner specified by the dynamics of the theoryJ23 and 
it is this response which is the topic of interest in the 
current paper. 

The experimental situation corresponds to applying an 
external magnetic field H , and then study the magnetic 
response of the superconductor to this field. Hence, a 
suitable thermodynamic description is coached in terms 
of a potential ^(H), which is a function of the inten- 
sive field variable H. In the simulations we have fixed n, 
which is analogous to fixing the magnetic induction, and 



a description based on the extensive field variable B is 
more appropriate. The tssp approaches are related by a 
Legendre transformationEZI. In principle the simulations 
could also be performed in an ensemble with fixed mag- 
netic field. Technically this would be achieved by adding 
the term 

H Lz2Tm/q 

to the action in Eq. |l|. This would promote n to a 
dynamical variable of the theory, and be more in ac- 
cordance with the experimental situation. However, a 
change n — > n ± 1 would require a global relaxation, and 
this would give very low acceptance rates, i.e. inefficient 
simulations. 

For type-I superconductors, superconductivity van- 
ishes for H > He- For type-II superconductors, a flux 
line lattice is formed at H — Hc^ , for smaller fields the 
magnetization in the sample vanishes due to the Meissner 
effect. By fixing n one can not study these effect directly, 
however it is possible to determine a corresponding field 
strength from n, see Ref. 

On the basis of simulations performed using the 
modified action Eq. |^, we have determined the effec- 
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tive temperature-renormalized interaction between two 
vortex-lines, and searched for the value of the GL pa- 
rameter, or more precisely its square, xi/u, where this 
interaction changes character from being effecti vely at- 
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tractive to being effectively repulsive. In section [I A 
have fixed n — 2 an d stu died the distance between the 
flux lines. In section 11 B we have generalized to real n, 
and used this to calculate a free energy difference be- 
tween states containing one and two vortices. This is 
also a measure of the sign of the effective vortex interac- 
tion, and hence an indication of whether we have type-I 
or type-II behavior. 

To obtain the results in Figs. |l|, and || we have 
performed simulations on cubic systems of size N — 
8,12,16,24,32,48, with /3g = 1.00. AU simulations 
have been performed in the broken symmetry state 
y < ijc, with particular emphasis on the values y = 
-0.04,-0.10,-0.20,-0.30,-0.40. For the two largest 
system sizes the final datapoints are averages of approx- 
imately 10^ sweeps, whereas approximately 10^ sweeps 
have been performed for the four smallest system sizes. 

The simulations leading to the results of Fig. ^ are 
quite different. They are performed for the fixed system 
parameters TV = 24 and f3c = 1-00, and for each value 
of m, we have performed from 2.5 • 10'' to 2.5 • 10^ MC 
sweeps through the lattice. One sweep through the lattice 
consists of (1) conventional local Metropolis updates of 
and A, and (2) gk^l radial updates of I'i/'l combined 



with overrelaxationr^ri of ip. 



r±{z) of these two vortices labelled by 1 and 2, and the 
average distance between them. 
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For type-I superconductors this distance should be in- 
dependent of system size, whereas for type-II we expect 
that this distance scales with the system size. Finite size 
scaling of d for various points in the (a;, y) phase diagram 
is shown in Fig. |l|. 
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A. Effective vortex interaction 

We first clarify what is meant by effective vortex in- 
teraction in this context. In the Ginzburg-Landau model 
at zero temperature, one may compute a pair-potential 
between two vortices which consists of an attractive part 
due to vortex-core overlap, and a repulsive part due to 
circulation of supercurrents (or magnetic fields) outside 
the vortex core. Ignoring fluctuation effects, this fur- 
nishes an adequate way of distinguishing between type-I 
and type-II behavior, by asking when the attractive core- 
contribution dominates the magnetic field contribution 
or vice versa. By effective interaction, we mean a ther- 
mally renormalized pair interaction which fully takes en- 
tropic contributions into account. At low temperatures 
the effective interaction will revert back to the standard 
pair-interaction described above, but will deviate as tem- 
perature is raised, and this is particularly relevant as the 
critical temperature is approached, as we shall see be- 
low. We also comment further on this in the Discussion 
section, where we elaborate on what we perceive to be a 
crossover between type-I and type-II behavior. 

In our simulations, the value of n has been fixed to 
n = 2. This corresponds to the case of two field- induced 
vortices which move around in the system under the in- 
fluence of their mutual effective interaction. During the 
simulations, we have measured the transverse position 



FIG. 1: The ensemble averaged separation (d) between the 
two vortices, as a function of system size. The two upper 
figures are indicative of type-I behavior, whereas the two lower 
ones indicate type-II behavior. 

In the part of phase diagram which we focus on, namely 
the region defined by the dotted line in Fig. 1 of Rcf . 0, the 
vortex lines are generally directed and almost straight, 
well defined line objects. This can be seen either by di- 
rectly taking snapshot pictures of the vortex-line con- 
figurations of the system, or by computing the mean- 
square fluctuations around a straight-line configuration, 
(|r^(z) — r^(O)P), for one vortex line. This is in contrast 
to the situation in the vicinity of the critical part of ydx) 
in Fig. 1 of Ref. |^, where the vortex lines loose their line 
tension via a vortex-loop blowoutE3E3. Consequently, we 
can consider the theory as an effective theory for an inter- 
acting pair of straight vortex lines which interact with the 
dimensionless potential V{d). If we make this assump- 
tion, the probability of finding the vortices separated by 
a distance d in a system of lattice points with periodic 
boundary conditions, is given by 



Pjv(d) 



e-^(^)l]jv(cl) 
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where J77v(d) is the number of configurations with a 
transverse vortex- vortex distance of d, and is just a 
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normalisation factor. riAr(d) can be calculated, either 
analytically in the continuum limit 




f-2arccos(|^)) 
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or by simple geometric counting in the case of a lattice. 
In the case of noninteracting vortices, i.e. V{<i) = 0, the 
expectation value of d is determined only by r2Ar(d), and 
we find the numerical value 



do = (d> 
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dd r2Ar(d)d w 0.38iV. 
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The separation do defined in Eq. [T^ will be used to es- 
tablish a numerical value of Xi^u. Namely, we can com- 
pute the averaged distance between vortices at fixed x 
varying y, or vice versa. In the latter case, we will use 
the criterion that if (d) exceeds som value cdo where c is 
some fraction, then we have type-II behavior, otherwise it 
is type-I. The quantity (d) at fixed y will turn out to be an 
S-shaped curve as a function of x, increasing from small 
values to large values as x is increased. We interpret this 
as yet another manifestation of the crossover from type- 
I to type-II behavior, and we have chosen to locate the 
crossover region Xi/u at the value of x where the curves 
change most rapidly, which is roughly when (d) w do/2. 
As we shall see (see Fig. |^), different crossover criteria 
give consistent results. The quantity Pwi^) can be es- 
timated from histograms, and then we can use Eq. || to 
determine the pair potential. Depending on whether we 
consider type-I or type-II superconductors we expect to 
see an attractive or a repulsive potential. Fig. || shows the 
potential V{<i) for the same points of the phase diagram 
as Fig. 



B. Free energy 

In Eq. I we have used n to indicate an integer number 
of flux tubes, but in principle there is no reason to limit 
n to integer values, and we will use S[a,ip,m] to denote 
a generalisation to real n. 

We have considered the free energy difference between 
a state containing zero vortices, i.e. m = and a state 
containing n vortices. We can not measure absolute val- 
ues of the free energy, but by differentiatingcZI 



-F(m) ^ T|.g-S(m) 



(11) 



with respect to m, and then integrating up to n, we can 
calculate AF{n) = F{n) - F{0), 
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FIG. 2: The effective interaction potential between vortices 
\^(d) as determined from Eq. |^. Observe the difference in 
vertical scale, in the lower panels (type-II) the interactions are 
much weaker than in the upper panels (type-I). The graphs 
correspond to the same points in {x,y) as those in Fig. |l| 
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FIG. 3: W{m) The straight line corresponds to 2m/NxNy 
which according to Eq. |l^ should be satisfied for m integer. 



To calculate A_F, we have then varied m in steps of 
Am — 0.05, and performed the integration in Eq. |u 
numerically. Using shift symmetries, it can be showrEI 
that W{n) is equal to 



W{n) 



2n 



(13) 



and the behavior for intermediate real values is shown 
in Fig. ||. Increasing n from to 1 costs a free en- 
ergy Ai^(l), and adding two vortices costs an amount 
Ai^(2). We wiU always have AF(2) > AF(1), but the 
question is whether AF(2) ^ 2Ai^(l). We may regard 
F(n-|-2)-|-F(n) — 2i^(?T,-|-l) as the discrete second deriva- 
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tive of the free energy with respect to particle number, 
which is nothing but the inverse compressibihty K^^ 
of the vortex-system. In the thermodynamic hmit this 
quantity can never become negative. However, its van- 
ishing signals the onset of phase- separation of the vortex 
system, which we again interpret as a lack of stability of 
the vortex-loop tangle, characteristic of type-I behavior. 



C. Vortex compressibility, separation, and 
crossover 

We next define a quantity AT by the relation 



AT = 



AT(1) Ai^(2) 
T, 2Tx 



(14) 



which means that AT measures the relative free-energy 
difference between adding one vortex to the system and 
half of that adding two vortices to the system. Intuitively 
it is therefore clear that it measures the sign of the vortex 
interactions, and hence determines whether we are in the 
type-I or type-II regime. AT > signals attractive in- 
teractions, i.e. type-I behavior, whereas AT < signals 



repulsive interactions, i.e. 
culated AT by using Eq. 



t^ 



pe-II behavior. We have cal- 
and Eq. n3, the results are 



shown in Fig. ^. The main qualitative result from these 
simulations is again that xi/ii{y) is a declining function 
of y. Note also that = —AT, and hence a posi- 

tive AT clearly implies phase-separation and instability 
of the vortex system, characteristic of type-I behavior. 
This is precisely what we see for small a; in Fig. ^ 
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FIG. 4: AT(a:) for different values of y. The attraction 
vanishes for xi/n < 0.5, and Xi/ii(y) is an increasing function 

of 

Finite size scaling of (d) and studies of AT differentiate 
nicely between strongly type-I and type-II superconduc- 
tors, but it is difficult to locate a value of xi/ii{y) with 
any great precision. Fig. || shows (d)(x) for different val- 
ues of y, along with a horizontal line at do/2, where do 
is the average separation between vortices had they been 



non-interacting. We have found that do ~ 0.38iV in our 
simulations. We have, rather arbitrarily, taken the inter- 
ception of this horizontal line with the curve (d)(x) as 



Vortex separation 




FIG. 5: The ensemble averaged distance between a 
pair of vortices, (d), as a function of the square of the 
Ginzburg-Landau parameter x — , for various values of the 
temperature-like variable y. The horizontal line is at do/2. 
Increasing y amounts to increasing the temperature. 

The curves of (d) [x) do not get significantly sharper 
with increasing system size, and there are no particu- 
lar sharp features in S[a,ip,2] as x is increased beyond 
Xi/n- Fig. ^ shows the intercepts from Fig. ||. Due to the 
features in the curves of Fig. ^, and how the results of 
Fig. ^ are obtained from them, we tentatively conclude 
that the computed line of Fig. ^, corresponding to the 
dashed line of Fig. 1 in Ref. is a crossover and not a 
phase transition. However, we comment further on this 
in the concluding section. 

As already indicated, there is some arbitrariness in the 
location of xi/ii{y) in Fig. |[ however the four points 
labelled by (lA, IB) and (HA, IIB) clearly are in the type- 
I and type-II regimes, respectivly. This is demonstrated 
in Figs. |l| and |. 



III. DISCUSSION 

From Figs. I and I we conclude that there is a crossover 
line separating effective attractive vortex interactions 
from effective repulsive ones, i.e. type-I and type-II. This 
line can either be crossed by changing x, i.e. lA — s- HA, or 
by changing the temperature i.e. IB — > IIB in Fig. ^. This 
means that for x values in a suitable range, we can have 
in principle have a temperature induced crossover from 
type-I to type-II superconductivity. Finally, we note that 
2^1/11 (y) deviates significantly from the mean-field value 
of xi/ii = 0.5. 

Deep in the type-I regime, we find clear evidence of at- 
tractive interactions. In the type-II regime the repulsive 
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FIG. 6: The computed crossover line yi/ii{x) separating 
type-I and type-II superconductivity. The black circular filled 
points are given by the intercepts between do /2 and the curves 
in Fig. H The four points labelled by (lA, IB, IIA, IIB) are 
the ones that were considered in detail in Figs. ^ and &_^The 
point marked T is the tricritical point in Fig . 1 of RefTM, see 
the discussion of finite size effects in section |III A| . The solid 
line connecting the points is purely a guide to the eye. The 
computed line above corresponds to the dotted line of Fig. 1 
in Ref. Qin the vicinity of (stri, J/tri)- 



interactions appear to be weak, and the results are es- 
sentially also consistent with two randomly placed vortex 
lines, i.e. not interacting, but not consistent with an at- 
tractive force between the vortices. Therefore, what the 
results unequivocally show is that by fixing the material 
parameter x and varying the temperature-like variable y, 
the character of the effective pair-potential is altered in- 
side the superconducting regime, significantly away from 
the critical line. 



A. Simulations 

The simulations with n — 2, and the simulations with 
a real me [0, 2] are quite different, and we will discuss 
them in turn. An important feature of all the simulations 
in the present work is slow dynamics. 

For the n = 2 simulation, where we have monitored 
d, we find that deep in the type-f regime the simulations 
are quite straightforward, the vortices stay close together 
with only small fluctuations, and a moderate number of 
sweeps is sufficient to get good statistics. However when 
we increase x towards xi/u the effect is not that d sta- 
bilises at a higher value, instead we get get fluctuations 
between a type-I like state where the two vortices are 
close together, and a type-II like state with large vortex- 
vortex separation. This picture persists as x is increased 
into the type-II regime, the only difference is that the 
fraction of time spent in the type-I like state state de- 
creases. In fact the results of Fig. ^ and Fig. || are in 
the type-II region quite close to what we would get from 



noninteracting vortices. What makes the simulations dif- 
ficult is that moving a vortex line in the transverse direc- 
tion is a global change, and thereby very slow. Timeseries 
of d show characteristic time scales of 10"* sweeps, so long 
simulations ^ 10^ sweeps over the lattice are required to 
obtain acceptable accuracy. A truly high-precision deter- 
mination of xi^ii(y) would surely benefit from a special- 
ized algorithm for the MC updates. 

To get good results one should take the — > oo limit. 
The conclusions from the results of Fig. |^ are based on 
this limit, whereas those drawn from Figs. |^,^ and ^ are 
based on the fixed system size N = 48. We have not 
performed a systematic study of finite size effects, but 
the curves in Fig. || do have subject to finite size effects 
in them. The trend is that curves move to the right 
upon increasing system size, this very likely explains the 
apparent discrepancy between the tricritical point (where 
the N oo limit has been applied), and the remaining 
points in Fig. ^. 

Note that Eq. ^ is a continuum field theory, and 
xi/ii{y) is not a critical point, hence the continuum limit 
/3g ~* oo should be taken. Our experience from the 
large-scale simulations performed in Ref. |^ indicates that 
(3g = f -00 provides conditions in the simulations already 
quite close to the continuum limit. We have therefore 
chosen to work with (3g = 1-00 and focused our efforts 
on considering large systems and long simulations. 

The Monte-Carlo computations of AT have been even 
more time consuming, because we have had to do the 
simulations for 41 different values of to. We have there- 
fore limited ourselves to considering only the system 
L = 24, Pg = 1-00; for a discussion of finite N j finite /3g 
effects see Ref. The relaxation time for these simu- 
lations has been particularly long in the limits to ^ f ^ 
and m — > 2~^ , and we therefore have performed much 
longer simulations in these limits than for intermediate 
values of to. From Fig. ||and Eq. ^it is seen that W(m) 
is a quantity of order 0[N^'^) whose finiteness originates 
in the difference between two 0(1) quantities. Conse- 
quently, it is difficult to get numerically precise results. 
This has in particular been the case in the limits to ^ 1~ 
and TO, 2^. 



B. Crossover 

The physical picture that emerges for the thermal 
renormalization of the vortex interactions, is the follow- 
ing. At low temperatures, for the K-values we consider, 
the system is in the type-I region with a fairly deep min- 
imum in the effective pair-potential between vortices at 
short distances, leading to an attractive interaction. At 
very short distances, we find a steric repulsion on the 
scale of the lattice constants in the system due to the 
large Coulomb barrier that must be overcome to occupy 
a link with two or more elementary vortex segments. This 
lenght scale represents the size of the vortex core in the 
problem. Upon increasing the temperature to the vicinity 
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of the line ydx), we do not find large transverse mean- 
derings of the individual vortex lines as we move along 
each vortex line, rather the vortex lines are essentially 
straight. Therefore, we believe that it is not entropic re- 
pulsion due to the bare steric repulsion in the problem, of 
the type which it seems reasonable to invoke for strongly 
fluctuating elastic stringsE2lEil that renormalizes the vor- 
tex interactions in the way that is seen in Fig.|. Rather, 
what appears to happen is that the vortex lines slosh 
back and forth in the minimum of the effective potential 
well as essentially straight lines. Hence, to a larger and 
larger extent as temperature is increased, they experience 
the hard wall in the interactions at small distances, and 
the weak attraction at large distances, effectively wash- 
ing out the minimum in the potential, thus making it 
effectively more repulsive. 

This is also seen in our simulations (not shown in any of 
the figures) when we monitor the transverse meandering 
fluctuations of each vortex line, {\r±{z) — r±{Q)\'^), as 
well as the mean square fluctuations of the intervortex 
distance, (d^) — (d)^, where d is defined in Eq. (0). The 
former is small deep in the superconducting regime, and 
remains small as the line ydx) in Fig. 1 of Ref. ^ is 
approached, while the latter increases dramatically as the 
dotted crossover line is crossed. It is precisely this fact 
which makes the simulations extremely time consuming. 

One should however keep in mind that, since we are 
considering the full GL theory in our simulations, and 
not the linearized London limit, it is in any case a drastic 
simplification to view the effective interaction between 
vortices as a simple pair potential. 

Finally, we note that, although our present simula- 
tions, which by necessity are on finite-sized systems, in- 
dicate that the change from type-I to type-II behavior is 
a crossover, we cannot rule out the possibility that it is 
elevated to a true phase-transition in the thermodynamic 
limit. More work is needed to clarify if this is indeed the 
case, but this will have to await the next generation of 



massive parallell computers. Questions that need to be 
addressed in this context, are: What is the order param- 
eter of such a transition, and what symmetry, if any, is 
being broken. 



IV. CONCLUSION 

We have considered the effective interaction between 
two vortices in the full GL model, and how this effec- 
tive interaction is influenced by thermal fluctuations. We 
have included fluctuations in the gauge fields, as well 
as the phase- and amplitude-fluctuations of the complex 
scalar matter field of the problem. We have found that 
the effective interaction changes from being attractive 
to being repulsive at Xj/n- This means a change from 
type-I to type-II behavior. We have found that xj/n is 
below the standard quoted value of 0.5, and is a func- 
tion of the temperature-like parameter y. This means 
that at the critical point, the value of the GL parame- 
ter that separates type-I from type-I behavior is smaller 
than I/V2. The line xi/ii{y) appears to be a crossover, 
and not a true phase transition. The above seems to 
offer a simple explanation for the experimental obser- 
vation that elemental Ta and Nb superconductors show 
a crossover from type-I to type-II behavior as the tem- 
perature is increased towards Tc- Previous explanations 
based on mean-field theories and not involving thermal 
fluctuations required two additional temperature depen- 
dent K-values to be definedO. 
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